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Abstract 

"vb" . 

Lattice QCD studies of hadron-hadron interactions are performed by computing the energy 

o ' 

{SJ ' levels of the system in a finite box. The shifts in energy levels proportional to inverse powers of 

the volume are related to scattering parameters in a model independent way. In addition, there 



are non-universal exponentially suppressed corrections that distort this relation. These terms are 
proportional to e~ m7lL and become relevant as the chiral limit is approached. In this paper we 
report on a one-loop chiral perturbation theory calculation of the leading exponential corrections 
in the case of / = 2 irir scattering near threshold. 
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I. INTRODUCTION 



Except in the case of infinitely heavy baryons, where an adiabatic potential can be defined, 
the interaction between two hadrons is studied with lattice QCD by numerically calculating 
energy levels of the system in a finite box. This is because in the infinite volume limit and 
away from kinematical thresholds, the two-hadron Euclidean correlator gives no information 
about the Minkowski space amplitude [1]. The alternative is to consider the system in a 
finite box, as is the case with numerical calculations anyway. The energy levels of a system 
composed of two hadrons are not simply the sum of the energies carried by each hadron, but 
there is an additional (usually small) shift that arises due to the interaction between them. 
The smaller the box, the larger the shift in energy levels. This volume dependence is inversely 
proportional to the volume and furthermore, there is a relation between the energy level shifts 
and the scattering phase shifts [2-4]. This relation, valid for energies below the first inelastic 
threshold is a consequence of unitarity and is thus model independent. 1 In addition to this 
power law shift in the energy levels, there are exponentially suppressed corrections which are 
not model independent and are the analogue of the exponentially suppressed corrections to 
the mass, decay constants, etc., in the single-hadron sector [5-9]. These exponential volume 
effects arise because the off-shell propagation of intermediate states is altered by the presence 
of the finite box, which allows them, for instance, to "wrap around" the lattice. As such, 
these effects are dominated by the lightest particle, the pion in QCD, and are proportional 
to e~ m ^ L with m-x the pion mass and L the linear dimension of the box. For simulations 
done with small enough quark masses such that the pions are within the chiral regime, these 
soft pion effects can be computed using the chiral perturbation theory (%PT) [13, 14]. The 
7T7T scattering phase shifts have been computed using lattice QCD following the universal 
finite volume method mentioned above [15-31]. As the chiral limit is approached [31] and 
more precise calculations appear, these exponentially suppressed corrections will need to be 
understood. 

Our goal in this paper is first to show the modification of the universal scattering formula 
for a hadron-hadron system in a box due to the exponentially suppressed finite volume 
corrections. Second, we compute the dominant exponential volume dependence explicitly 

1 By model independent relation we mean a relation valid whether one is considering QCD or some other 
theory, as long as this theory obeys unitarity, locality, etc. 
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for the two-pion system in I — 2 channel near threshold by use of the leading loop-order 
two-flavor xPT. 

II. FINITE VOLUME vrvr SCATTERING 

A. Power law and exponential volume dependence 

As discussed above there are two types of volume dependence of the energy levels of two 
hadrons in a box: power law (proportional to 1/L 3 ) and exponential (proportional to e^ m ^ L \ 
The first is exploited by the finite volume method to extract information about scattering 
parameters [2-4]. The second, usually numerically smaller, appears as a correction to the 
relation between energy levels in a box and scattering parameters. In order to compute 
the exponentially suppressed terms we need to separate them from the larger power law 
contribution. 

Figure 1 shows all ttit scattering diagrams which contribute at one-loop order. As we will 
discuss in more detail in the next section, the power law corrections arise only from s-channel 
diagrams as shown in Fig. 1 (a), where the intermediate particles can be on-shell, and thus 
propagate far and "feel" the finiteness of the box. In all other diagrams the intermediate 
particles are very off-shell, cannot propagate farther than a distance of order l/m w and 
therefore have only small, exponentially suppressed sensitivity to the size of the box. 

Let us now discuss the general form of finite volume corrections. Consider first the 
pion propagator at finite volume which is a function of the spatial momentum k = 2nn/L 
and the energy E. It will have poles for values of E corresponding to the values of the 

— # 

energy of a pion in the box. In particular, for k = the pole will be at m„(L) (the "finite 
volume mass"), differing from the (infinite volume) mass of the pion by an exponentially 
small quantity proportional to (m 2 / (47r/) 2 ) e~ m7rL / LV mL [6, 7, 9]. The extra suppression 
factor m^/(47r/ 7r ) 2 is due to the fact that only loop diagrams contribute to the finite volume 
corrections. 

The volume corrections for systems with more than one hadron are more subtle. The 
reason is that there are two kinds of volume corrections to the energy levels: a power law 
one described by the Luscher formula and the exponentially suppressed ones. To understand 
how to separate them let us first look at the infinite- volume, S-wave scattering amplitude, 
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T(s), with on-shell external pions. It is given at one loop by 2 : 

T(s) ~ T(°)( S ) + T«( a ) + T«( S ) + ir#( S ) 

(TW( S )) 2 



(1) 



r(^( a )-r«( a )-r«( a )-£ri5( s 

where T^ n ^ is n-th loop contribution, the s-channel contribution is separated into its real 
part, T^(s), and imaginary part, T s j (s), and all other contributions at one-loop including 
t- and w-channels are denoted by T^(s). The imaginary part, T s j(s), comes from picking 
in the loop integration, the particle poles in both pion propagators, such that both the loop 
pions are on-shell. The loop integral is then proportional to the phase space volume and is 
given by 
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The fact that the imaginary part is determined by the tree level amplitude is a consequence 
of the optical theorem. 






(a) (b) (c) 





(d) (e) 

FIG. 1: The one-loop diagrams which contributing to the tttt scattering amplitude. Only the s- 
channel diagram, (a), contributes to the power-law volume dependence. Diagrams (b) and (c) are 
the t-, and u-channel diagrams, respectively, while diagram (e) represents wavefunction renormal- 
ization. All these diagrams contribute to the exponential volume dependence. 



2 We are considering the s-wave projected amplitude and disregarding the mixing with higher partial waves 
induced by the breaking of rotational symmetry. 
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It is useful to define the i^-matrix [10], which at one-loop is given by: 



Since the scattering amplitude can be written in terms of the phase shift, S(s), as 

= 327r^ 1 = 32^ e^)-l 

y/s-4ml cotS(s)-i y/s - 4ml 2i ' 
the relation between the X-matrix and the phase shift is then given by 



1 ' I"' lm ^cot5(s). (5) 



K(s) 32tt V s 

Now we look at the finite volume amplitude T(s) 3 . It is computed in the same way as 
the infinite volume amplitude, except the loop integrals are substituted by sums over the 
momenta allowed in a finite box. The important point to keep in mind is that sums where 
the summand is regular are, at large enough L, well approximated by the analogous integral, 
up to exponentially small terms. If the summand, however, contains a singularity, power law 
dependence on the volume arises. As mentioned before, only the kinematics of the s-channel 
diagram allows for both of the intermediate pions to be on-shell simultaneously. This implies 
that the summand in the sum over the loop momentum contains a singularity and leads to 
power law volume corrections. For the remaining diagrams no singularities are present and 
only exponentially suppressed corrections can arise. As it will be shown explicitly below, the 
finite volume amplitude then has the following form: the tree term remains the same, the t-, 
u- and the real part of the s-channels pick only exponential corrections but the imaginary 
part turns into the term with power law L-dependence (and is real at finite L): 

T(s) * r<«) (s) + Tfi'W + 7f» + AT«( S ) + g£) 

r(o) (s) _ r« (s) - r« (s) - at«( s ) - s (igf ) ' 

where s = 4(k 2 + m%), AT e xp(s) is the finite volume correction to T^(s) +T^(s), and S is 
a universal (independent of the interaction) function of s [11, 12]: 

5 O - 4 ^ ^ £ - / (&1 ^ = tL E ^ - ^ m 

- ^ V 7 |n|<A n 4tt 2 



3 By finite volume scattering amplitude we mean the amputated four-point correlator since, of course, there 
is no scattering at finite volume. 
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Combining Eqs. (1,3,5,6) we have 

T(s) ~ _ j; 



AT^l(s) i c( {s-Aml)L^ \ 

K(s) (T(°)(s)) 2 167^1,^ ^ 16tt 2 ^ 

16^^ 

fccot<5( S ) - I6n V~s0f^ - (^g^) 



(8) 



The energy of the states in the box are determined by the location of the poles of the finite 
volume amplitude determined by the solution of 



(T(°)(s)) 2 ttL V 16vr 2 
We recognize in Eq. (9) the familiar form of the Liischer relation modified by the finite volume 
correction: the quantity — 1 6 7T EW(sf) 2 * s ^ ne sou ght-after (exponentially small) correction 
to k cot S(s). In this work, we will focus on the 2-pion correlator near threshold in the center 
of mass frame, for which the infinite volume energy is given by ^/s = Im^. The solution of 
Eq. (9), s*, will be away from threshold by an amount given by y r s* — 2m 7T w y r s* — 2m 7r (L) ~ 
1/ f 2 L 3 . Therefore, for s pa s* the correction term, A(/ccot5(s)) = — 1 Qit^fs ^0)^2 , can be 
approximated by — 327rm 7r ^7(4^?^ , the difference being suppressed by ~ 1/L 3 . 

It is customary to expand Eq. (9) in powers of k 2 ~ 1/L 3 . Up to the first three orders 
of this expansion (near threshold), A;cot<5 can be approximated by the inverse scattering 
length, 1/a, resulting in 

^- 2m ' = ^( 1 + Cl z + < Hz) +'■■)• (10) 

where c\p are known numerical factors. A generalization of this formula including the expo- 
nentially suppressed corrections, however, is not useful. The error in using Eq. (10) instead 
of Eq. (9), that is, the error in the extrapolation from s* to 4m 2 , is of order 1/L 3 , which 
is parametrically larger than the exponential corrections we are interested in. Numerically, 
it may be the case that, for a set of simulation parameters, the exponential term is larger 
than the 1/L 3 terms. But the analogue of Eq. (10) one would obtain by formally counting, 
for instance e~ m7rL ~ 1/L 2 , would involve the effective range in addition to the scattering 
length. In any case, it is unclear a priori that a set of simulation parameters exist where 
this kind of expansion is useful and we will not pursue this line of thought in this paper. 
We now will compute the exponential corrections to k cot S(s). 
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B. The 7T7T scattering amplitude 



The 7T7T finite volume correlator in the 1 = 2 isospin channel, for arbitrary momentum 
(in the chiral regime) is given by 

(3s -2m 2 - £li p?) 



% = — 



f 2 



+ 



10s 
9P 



6m^ 



4 



+ 
+ 
+ 
+ 
+ 



4 2/ \ v f \ v ~\ 

^pPtPt + gp[P t + 3 (^ + P^J [P* + 3 & 2 + P*>) \ WviPt) 

4 2/ \ ^ / \ "1 

L g^P«P« + 777i ( P« + 3(pi + Pa) ) (p« + 3(p 2 + Ps) ) iJ^(p u ) 

Ms- 



9f 2 

3m 2 ) 2 . — 4 



9/ 2 

4-um 2 



m 



Atm 2 



+ 



2t 2 



m 

7 2 



2v*_ 

3f 2 + 3f 2 



4(m 2 _-w) 



iJ(p u ) - 



3f 2 3/ 2 . 
8(s-3m 2 ) 



iJ(pt) 



4K - 1) 



9/ 2 



P^MPs) + 3/2 ' rfiJM 



PuiMPu) ~ [(t 



2i2 

p 



2m 2 ) 2 + (u 



2(s - 2m 2 ) 2 + (f - 2m 2 ) 2 + (u - 2m 2 ) 2 



2\2 



2m 2 ) 



32 m 4 

T 'T 2 " 



(ii) 



In the above expression, m is the (volume independent) tree level pion mass, / 132 MeV 
is the (volume independent) tree level decay constant and the £j's are the Gasser-Leutwyler 
coefficients of counter terms appearing in the chiral Lagrangian at next-to-leading order 
(NLO) [13]. The loop integrals/sums are given by J, J(P), J^P), and J^ V (P) and will be 
defined below shortly. The external momenta p i: % — 1 • • -4 are described in Fig. 1, and the 
Mandelstam variables are employed: s = p 2 with p s = pi + p 2 , t = p 2 with p t = pi — p 3 , 
and u = p 2 with p u = p\ — p 4 . In the above equation, the first term is the leading-order 
(LO) tree level contribution, and the remaining terms come from the one-loop diagrams 
shown in Fig. 1 and tree level diagrams with vertices of the 0(p 4 ) Lagrangian [13, 14]. In 
the NLO contributions, we have approximated p 2 = m 2 , as the corrections to this at finite 
volume are beyond the order we are working. In the one-loop terms, we have expressed 
all contributions in terms of the bare pion mass as the difference is higher order in the 
chiral expansion. One can, of course, choose to express the scattering amplitude in terms 
of the "lattice quantities" such as m n (L) and f n (L), which are measured directly from 
Euclidean correlation functions by lattice simulations. Converting the bare quantities into 
lattice quantities involves additional tadpole loops, which will affect the form of A%!xp , the 
finite volume corrections to the two-pion amplitude. The exponential volume dependence 
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of course doesn't depend upon whether one expresses the amplitude in terms of either the 
bare of physical parameters, and so it is useful to use the form which is simplest, that in 
terms of bare parameters. In what follows, we will be interested in the 2-pion correlator 
near threshold, for which the external pion momentum are given by Pi ~ {\\f~s, 0). 



C. Loop integrals/sums at one-loop 

The loop integrals/sums appearing in Eq. (11) are defined by 

[ <!S±L V 1 (12) 



J 



j(p) = [ v - 1 as) 

Uy } J 2ttL 3 ^g 2 -m 2 (P + g) 2 -m 2 ' 1 ' 

J(P)= - (14) 

^ M ' J 2ttL 3 ^ g 2 -m 2 (P + g) 2 -m 2 ' 1 ' 

L 



and 



J^(P) [ ^ ^3 q 2^ m 2 (p + q) 2 _ m 2 ■ ( 15 ) 



Note that an integral is taken along the th component whereas sums over discrete momenta 
are taken with cubic symmetry. Finite volume effects in the loop integrals/sums in Eqs. (12)- 
(15) can be computed by first evaluating the qo contour integral and then using the Poisson 
resummation formula, 

£ Sw? 5 ^" 4 "- (16) 

The difference between the finite volume and infinite volume loop integrals/sums is given 
by the second term in the right-hand side of Eq. (16), and is always ultraviolet finite. If 
the function f(q) is regular, this difference is exponentially suppressed in the large L limit. 
Power law dependence on L can however appear if f(q) has a singularity, i.e., the case when 

P = Ps- 

Let us now evaluate the difference between the finite and infinite volume integrals/sums 
given in Eqs. (12)-(15). We shall define the difference between the finite and infinite volume 
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integrals/sums as 



A/ S /(FV) -/«»)=/£ i.E"/ 



d 3 g 

L ■• ^ ./ (2tt) 3 

L 



(17) 



where it is implicit that we regulate the sum and integral in the same manner such that the 
UV divergences cancel. The tadpole integral in Eq. (12), which contributes to m n and /„. 
and the loop diagrams of Fig. 1, has the following volume correction; 



i 



r,3 ^2 



L 3 ^ J ( 2 7t) 3 

q ~ L 



d 3 q 



q 2 — m 2 



L3^~ J (27 



2u q ' 



m 



4tt 2 L ^— ' \n 



V ^KMmL). 



(18) 



where u q = \J q 2 + m 2 . The mass and decay constant measured in lattice simulations are 
thus given to NLO by [7] 4 

iAX' 



m\{V) = ml 
UL) = u 



1 + 



f 2 
2iAl~ 



m 



f 



_ iX(L = oo) + iAX M z m 2 
1 + — =r + 



P 



P 



. 2iX{L = oo) + 2iAl 2£ 4 m 2 



(19) 
(20) 



P P 

Using the asymptotic form of the Bessel function, one can see that for large ml, the volume 
shift of the pion mass is exponential [6, 7], 



Ami iAX 



ml f 2 25/%3/2 L/ 2 



m^ x - 
7\f2 2^ 



-n m^L 



c(n) 



3 / 2 Vm~j: 



n 



1 + 



8 nm^L 



+ ... 



(21) 



n=|n|^0 

where the ellipses denote more terms in the asymptotic expansion of the Bessel function 
and c(n) is the multiplicity factor counting the number of times n — \n\ appears in the 
3-dimensional sum. Note, this sum is not over integers, but rather over the square-roots of 
integers. In Table I, we list the first few values of the multiplicity factors. 

Power law L-dependence can only occur through the integrals/sums in Eqs. (13)-(15) 
when P 2 > 0. For the center-of-mass scattering kinematics we are considering here this can 



4 This relation is known up to two loops [9, 32]. 
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TABLE I: Here we list the first few multiplicity factors which arise when converting the three- 
dimensional sum to a scalar sum. 



n 


1 y/2 V3 v 7 ! y/5 


c(n) 


6 12 8 6 24 



only occur for P = p s , since p s = s > 0. As argued above we will only need the amplitude 
at threshold, i.e., p s = (2m 7r ,0) and p t = p u = 0, except for the terms with power law 
L-dependence. Consequently, we will need only the values of AJ(P = 0), AJ (P = 0), and 
AJ 00 (P = 0) for t- and w-channels as well as AJ(P = p s ), AJ (P = p s ), and AJ" 00 (P = p s ) 
for s-channels. The J integrals/sums at P = can be shown to be related to 1, giving the 
volume difference: 



iAJ(0) = — 



iAJ (0) = 0, 



d 3 q 
(2tt) 2 



1 



d 



dm 2 



(iAl) , 



(22) 
(23) 



and 



^AJ 00 (0) 



-y- 



L 3 



d?q 
(2^f 



+ 



dqo 
2tt 



iAl + m 2 iAJ(0) + 3 — iAl - -m 

\ 6 3 dm 1 

-iAl. 



m 



+ 



q 2 — m? (q 2 — m 2 ) 2 (q 2 — m 2 ) 2 



2,A^ , of d /AJ 



(24) 



The power law volume dependence appears in the remaining integrals/sums. In those we 



keep s away from the threshold value and take ^/s = 2\/k 2 + m 2 . After performing the qo 
integral, we separate the singular piece of the summand from the rest as 



iJiPs) 



1^11 l^i l ^ 



4L 3 ^ uj„q- 
q 



4u k L 3 ^ q- 



1 Uq~ UJ k 

k 2 4L 3 cJgCJfc g • 



A' 2 



(25) 



The first term contains a singularity when the internal momentum coincides with the external 
momentum, while the second term is regular. The difference AJ(p s ) is then 



iAJ{p s ) 



S 



k 2 L 2 



n^0 



iLq-n Uq-UJ k 1 

<f 2 - A; 2 4w fc cj„ 



(26) 



V 

iAJ cxp (p a ) 
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The first piece above is the promised universal function containing the power law volume 
dependence. The summand in the second term contains only exponential finite volume 
corrections. This term, contributing to A%^, can be computed at the s = 4m 2 , k = 
threshold point, 

1 1 ^ 1 I" 00 , ylme i2ny W 



iAJ exp (p s ) 



E^ 



167T 2 Ly/m 2 + k 2 tri I'" 



dy 



y 2 + ^ 



y 2 + ^ + 



k 2 L 2 m 2 L 2 



167T 



L E 



n^O 



i^od^l mL)L_i{\n\ mL) + mL)L (\n\ mL) — — 



4vr 2 
1 



4vr 2 



n\ mL 



(27) 



where L v is the Struve function. To get the second line of Eq.(27), we have neglected terms 
which are suppressed by 0(k 2 /m 2 ) relative to the first. For the two-pion system, this is 
approximately given by ^ ~ <C 1. The asymptotic expansion of iAJ exp (p s ) is given 
by 



iAJ cxp (p s ) ~ 



(4tt) 2 (ml) 3 / 2 



E 



-n m w L 



c[n) 



n 



3/2 



1 - 



8nm n L 



+ 



(28) 



n=|ra|^0 

Again, one can see that these volume corrections to the integral are exponentially suppressed. 

The finite volume dependence of the other s-channel loop integral functions, iA,Jo(p s ) 
and iAJqq{ps) become simpler to evaluate by first observing that the summands can be 
separated into the following pieces: 



QQ-^n (90 +Pso) 2 -U 2 



1 + 



1 



(Psof 



1 



g 2 



L0„ 



(PsO + %) 2 - 
1 



il-ul (qo+Pso) 2 -u 2 



ll-^V (% + Pso) 2 ~ 



One then obtains, 



and 



iAJ (p s ) = -^AJ( Ps ) 



(29) 



iAJ 00 (p s ) =iAI 



l_ yl 



1 + 



k 2 



1 s 

-iAX + - iA J(p s 



(30) 



Having these tensor integrals/sums written in terms of the scalar integrals/sums iAJ{p s ) 
and iAX and using the scattering amplitude in Eqs. (1) and (11), one can now verify that 
the coefficient of «S(^r) in the amplitude is what was promised in Eq. (1). 
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At' 1 ' (4m 2 ) 

FIG. 2: Ratio of the magnitude of the exponential correction term A(k cot 5) = —32irm n ^(o^m 2 ^2 
to A; cot 5 for different values of the pion mass. 

D. Exponential volume correction to the I = 2 -kit correlator 

Collecting the results for the sums/integrals in Eqs. (18-30) and using the amplitude in 
Eq. (11) we can now compute the correction term in Eq. (9) for 1 = 2 two-pion system near 
threshold. We find 



A(JfeC0t<5(s)) = -327T7B 



ATexp(4m^) 



57T 



"(T(°)(4m2))2 m n [3 71 ' 3m\ 

227 1 



^-iAI + ml^-iAl + 2iAJ cxp (4ml 



\ „ , e nm * L 



c(n)- 



^2tt 7^, n ' ' y/nm^L 

n=\n\ftO 



24 nrriirL 



(31) 



Equation (31) is our main result (the first line being the exact one-loop answer and the 
second line the asymptotic expansion in m n L). In this expression, one can use either the 
bare parameters, the physical parameters or the finite volume parameters as the difference is 
higher order than we work in either the chiral expansion, or in the exponential dependence. 
It is most convenient to use the values oim n (L) and f n (L) directly measured in a given lattice 
simulation. In Fig. 2 we plot the ratio of A(k cot S(s)) to the one-loop value of k cot S(s) 
using Eq. (31) as a function of L for some reasonable values of m n . We find the finite volume 
corrections to be relatively small, a few times smaller than the statistical and systematic 
errors quoted in recent simulations. An error of about 10% was quoted in reference [31] 
for the determination of the scattering length for a pion mass of m n ~ 290 MeV and a 
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box size of L ~ 2.5 fm. 5 The finite volume correction from Eq. (31) for these parameters 
is approximately 1%. These corrections however grow fast with the approach to the chiral 
limit, and they become non- negligible as smaller pion masses are used and statistical errors 
are reduced in simulations. 

III. DISCUSSION 

We have described the leading exponential volume dependence expicitly for the scattering 
parameter of a / = 2 two-pion system near threshold in a box, by extending the one-loop 
%PT calculation of pion scattering [13] to include the volume dependence. The exponentially 
suppressed volume corrections can distort the universal relation between the infinite volume 
scattering parameters and the power-law volume dependence of the two-particle system, 
especially as the chiral limit is approached. An important point we want to stress is that 
the useful way to add the exponential volume dependence to the relation between infinite 
volume scattering parameters and the energy of the two-particle system in a box, is via 
Eq. (9), which allows an understanding of the leading exponential volume dependence to 
k cot S(s). This is contrast to the notion of studying the exponential volume dependence of 
the scattering length, the effective range etc., separately. 

It is important to stress the limits of validity of the present calculation. On one hand, the 
pion masses should be small enough so that the chiral expansion is converging. From the 
experience acquired in the three flavor case, where kaon loops are a borderline case for the 
convergence of the expansion, one expects chiral perturbation theory to be useful for m n < 

5 In Rcf. [31], Beane et. al. determined the I = 2 mr scattering length for various pion masses using a 
mixed action simulation with Domain- Wall valence quarks and staggered sea quarks [33]. Because of the 
mixed-action, the mesons composed of sea quarks and the mesons composed of valence quarks receive 
different mass shifts from the finite lattice spacing. This means that even when the sea and valence quark 
masses are tuned equal, there are still partial quenching effects in the simulation. In Ref. [34] , the partial 
quenching and lattice spacing corrections to the I = 2 mr scattering length were worked out for this 
mixed action theory. It was shown that the these two lattice artifacts were largely suppressed, and almost 
non-existent for the mass tunings used in Ref. [31]. However, as shown in Ref. [34], there are still partial 
quenching effects and in particular, in the t- and w-channel diagrams the hairpin contributions can be 
significantly more sensitive to the boundary effects. For / = 2 these effects are only exponential, and for 
the pion masses and box sizes used in Ref. [31], we have found they are the same order of magnitude as 
the corrections of this paper, and thus not-significant to the work of Beane et. al. These effects are being 
worked out in detail in Ref. [35] 
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500 MeV (of course the exponential volume dependence for a 500 MeV pion, or kaon will be 
negligible). Also, the box size has to be large enough so the usual power counting used here 
(the so-called "p-counting") is valid. When L is much smaller than the inverse pion mass, 
another power counting is required such as the e- [36] or the e'-regime [37]. Additionally, 
we have neglected corrections which occur from higher loops, all of which are suppressed by 
additional factors of (m 7r /47r/ 7r ) 2 and some of which are suppressed by additional exponential 
factors of e~ milL . The diagrams with this extra exponential suppression result from two- 
loop diagrams where intermediate states in both loops are purely off-shell and hence "going 
around the box" . 

We have focussed on the exponential corrections to phase shifts close to threshold. One 
can easily extend this work to include the exponential volume dependence of the phase 
shifts at higher energies. Alternatively, one can access non zero momenta by using twisted 
or partially twisted boundary conditions to probe the low-momentum dependence of the 
scattering amplitude [38-45]. This method may boost the entire two-particle system, how- 
ever, requiring the extraction of scattering parameters in a boosted frame [46-48]. Our 
methods generalize trivially to this case also. These methods can also be extended to other 
interesting two-hadron sytems [11, 12], where these exponential volume effects may be more 
significant. 
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